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1 Introduction 



Working out propagators is the difficult part about formulating quantum 
field theoretic perturbation theory on exotic backgrounds. It is typically 
accomplished by solving the differential equation that the propagator must 
obey, however, this procedure is ambiguous up to a homogeneous solution. 
It has long been realized that some choices for this homogeneous solution 
do not make the resulting Green's function into a true propagator. That 
is, the Green's function does not correspond to the expectation value of the 
time-ordered product of two free fields in the presence of any state pp. 

There is nothing mysterious about the problem, nor does it require field 
theory to understand. Consider the simple harmonic oscillator whose position 
as a function of time is q(t) and whose Lagrangian is, 

t 1 -1 1 2 2 /i \ 

L = -mq muj q . (I) 
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The propagator equation for this system is, 



— m 



(■^y + oo 2 ]tA(t;t') = i6(t-t') . (2) 



The general solution to this equation which is symmetric under interchange 
of t and t' has three free parameters, 

% r i 

iA(t;t') = sin u\t — t'\ + a cosiut) cos(ut') 

+/3sin u(t+t') + 7 sin(ujt) sm(ut') . (3) 

Although any choice of a, [3 and 7 in gives a solution to the propagator 
equation d2J), the result is not a propagator unless they obey two inequalities, 

1 X 

a + 7 > and aj > -3 2 . (4) 
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To see this, note first that the Heisenberg picture operator q(t) can be ex- 
pressed in terms of its initial position and momentum as, 

q(t) — % cos(ujt) + sm(ujt) . (5) 
muj 



For iA(t; t') to be a propagator there must be a state 
iA(t;t') = UT\q(t)q(t')} ^ 



such that, 



sin 



u\t-t'\ 



qoPo+poqo 



Imuj 



+ (ij) ij}^ cos(atf) cos(ut') 

ip) sin(cjt) sin(u;£') 



2m 2 u 2 



V>) sin cj(t+t') 
We can therefore identify the constants a, /3 and 7 as, 

9oPo+Po?o 



2mc<j 



For the ground state one has, 
Ground State = 



a 



■0 
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2m 2 w 2 
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and /3 = . 



(6) 



(7) 



(8) 



(9) 
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For a general state the first inequality in (jlj) results from requiring the ex- 
pectation value of the energy to be greater than or equal to the second 
is just the Schwarz inequality. 

A more subtle set of issues can arise in gauge theories. To understand 
the one of interest for this work we must digress to explain the difference 
between an "exact gauge" and an "average gauge" [2] . The former is obtained 
by choosing the gauge parameter to make the vector potential obey some 
equation at each point in space and time. This is the normal type of gauge 
fixing in classical field theory. Familiar examples are, 



Temporal Gauge 
Coulomb Gauge 
Lorentz Gauge 



A (t,x) = , 
V-A(t,x) = , 

d^A^t,x) = -A (t,x)+V-A(t,x) 



. 



(10) 

(11) 
(12) 



Although exact gauges can be used in quantum field theory the more common 
type of gauge fixing is accomplished by adding some noninvariant term to 
the invariant Lagrangian. For example, the Feynman gauge Lagrangian is, 



C = ~F lv F^-h^A ll ) 



(13) 



The functional integral representation for this type of gauge condition can 
be viewed as a weighted average of exact gauges. For example, the Feynman 
gauge functional formalism results from imposing the exact gauge, 



d*AJt,Z) = f(t,£), 



(14) 
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where f(x) is a (D-number field. One then functionally averages over f(x) 
with a Gaussian weighting functional, 



[df] exp[- l -Jdt Jd 3 xf 2 (t,x)] . (15) 



From this discussion it is obvious that a fairly involved set of functional 
changes of variables connects the exact gauge conditions of the canonical for- 
malism to the average gauge conditions typically employed in the functional 
formalism. The late Sidney Coleman worked this out explicitly for flat space 
on the manifold R 4 to derive the Faddeev- Popov ansatz for this case [3], 
but the result is often assumed without justification for general metrics on 
any manifold. We suspect that the unjustified use of average gauge fixing is 
behind a dispute about the graviton propagator on de Sitter background. 

It is easy to see that certain average gauges can be problematic when 
linearization instabilities are present. Consider flat space electrodynamics on 
the manifold T 3 x R. Because the spatial sections are compact, both sides 
of the spatially averaged, \i = Maxwell equation must vanish separately, 

O u F^ = J^ Jd 3 xdiF i0 {t,x) = Jd 3 xJ°{t,x) . (16) 

Because this zero charge constraint follows from the invariant field equations 
it must be true as well in every exact gauge. However, naively imposing an 
average gauge can result in a very different theory. For example the field 
equations of Feynman gauge (IT3"|) are, 

-d 2 t + V 2 ] A»{t, x) = J"(t, x) . (17) 

These equations can be solved for any total charge. One can argue about 
how the problem happened, or how significant it is, but there cannot be any 
doubt that something went wrong. 

The issues we have been discussing are relevant to a debate between 
cosmologists and relativists concerning perturbative quantum gravity on de 
Sitter background. From the perspective of inflationary cosmology it is natu- 
ral to view de Sitter as a special case of homogeneous, isotropic and spatially 
flat geometries whose invariant element in co-moving coordinates takes the 
general form, 

ds 2 = -dt 2 + a 2 (t)dx ■ dx . (18) 
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One gets the open coordinate sub-manifold of de Sitter by setting the scale 
factor to a(t) = e Ht with constant H. For any scale factor the transverse- 
traceless components of the graviton field obey the same equation as the 
massless, minimally coupled scalar [1], 



The power spectrum for gravitational radiation [5] is proportional to the 

3 

canonically normalized, super-horizon scalar mode functions u(t, k) ~ H/k?, 



From scale invariance — which would be exact in de Sitter — one sees that 
the mode functions of physical gravitons diverge too strongly at small k to 
give a convergent result for the Fourier mode sum of a part of the graviton 
propagator which must be present in any gauge. It follows that there can 
be no de Sitter invariant graviton propagator, just as there is no de Sitter 
invariant propagator for the massless, minimally coupled scalar [6]. 



People who abhor the breaking of de Sitter invariance typically dismiss 
it as unphysical, but this argument cannot be accepted because the tensor 
power spectrum is certainly physical. (Indeed, strenuous efforts are underway 
to observe it!) Nor is there any support to be gained from the tiny distinction 
between de Sitter and primordial inflation, which typically makes the infrared 
behavior worse in any case. So one would think that the noninvariance of 
free gravitons on de Sitter must be accepted as universally as that of the 
massless, minimally coupled scalar. This is is not so because relativists have 
been able to find de Sitter invariant solutions to the propagator equation in 
average gauges [Tj. Explicit solutions in what seem to be valid gauges are 
just as compelling as inferences from the tensor power spectrum. However, 
we have just seen that average gauges may not be reliable on manifolds such 
as de Sitter which possess linearization instabilities. 

The early recognition of a problem [8] with one of these de Sitter invariant 
solutions led to the development of a noninvariant, average gauge [HI HQ] . The 
associated propagator has been shown to obey the Ward identities at tree 
order [TTJ and one loop [12]; and the only fully renormalized, dimensionally 
regulated loop results for quantum gravity on de Sitter background have 
been obtained using it (131 [HI [15] . Although the gauge fixing functional of 
this propagator is not de Sitter invariant, it does preserve the one-parameter 




(19) 



V h ~Gx \u(t,k)\ 2 x k 3 ~ GH 2 . 



(20) 
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subgroup of dilatations, 

t — > t-^\n(k), (21) 
x — ► kx , (22) 

so the fact that the propagator breaks dilatation invariance cannot be blamed 
on the gauge. Moreover, Kleppe has shown that the propagator's de Sitter 
breaking is physical by the standard technique of supplementing naive de 
Sitter transformations with a compensating gauge transformation to restore 
the noninvariant gauge condition [16J. 

As with the cosmological power spectrum, one would think these results 
decisive, but the interest in a de Sitter invariant graviton propagator persists 
|17j . What seems to be necessary to settle the issue is two things: 

• A proof that the average gauges for which de Sitter invariant solutions 
have been found are not valid; and 

• An explicit construction of the graviton propagator in an exact, de 
Sitter invariant gauge which is valid over the full de Sitter manifold. 

Of course the imposition of a de Sitter invariant gauge would make the prop- 
agator equation de Sitter invariant, but that does not imply a de Sitter 
invariant solution for the graviton propagator any more than it does for 
the massless, minimally coupled scalar propagator which obeys an invariant 
equation but is not invariant [6]. If the graviton propagator in an exact gauge 
breaks de Sitter invariance then the breaking must be accepted as physical. 
This would not only resolve a contentious dispute, the resulting propaga- 
tor might be easier to use and it would reduce the number of counterterms 

OH eng. 

In section 2 of this paper we give a proof that the average, de Sitter in- 
variant gauges are not valid; in subsequent sections we develop the machinery 
for constructing the propagator in exact, de Donder gauge. The technique 
for our construction is to find the field-dependent gauge transformation ^[h] 
which enforces the gauge condition, 

h'^ = V - 2£ K *[h] such that 9 pa [h', p , a - \h' p ^] = . (23) 

(In this and subsequent formulae, g pu stands for the spacelike background de 
Sitter metric and a semi-colon denotes covariant differentiation with respect 
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to this background.) Then we use this transformation to carry the existing, 
noninvariant propagator [HJ HO] into de Donder gauge. If the de Sitter break- 
ing of the existing propagator is completely due to the noninvariant gauge 
then the de Donder gauge result should be invariant; if the de Sitter breaking 
is physical then it should persist after the transformation. 

As a warmup for de Sitter gravity we here carry out the same exercise 
for de Sitter electromagnetism. That is, we find the field- dependent gauge 
transformation 9[A] which imposes exact Lorentz gauge, 

A'p = A^ — d„9[A\ such that d^(^g^A' v ) = . (24) 

Then we use this transformation on the photon propagator in a noninvari- 
ant, average gauge [TQl HE] • Because we already know the unique, de Sitter 
invariant solution for the propagator equation in Lorentz gauge [19] , obtain- 
ing this known solution by transformation demonstrates that the technique 
will remove de Sitter breaking that arises from using a noninvariant gauge 
condition. The simpler setting of electromagnetism, and the close relation 
between the noninvariant electromagnetic and gravitational gauges should 
also teach us much of value for the main project. 

This paper consists of six sections of which the first is ending. In sec- 
tion 2 we review the functional changes of variables which carry one from an 
exact, canonical gauge to a covariant, average gauge, with special attention 
to what goes wrong when linearization instabilities are present. The con- 
text is flat space electrodynamics on the .D- dimensional manifolds R D and 
T D ^ 1 xR. In section 3 we switch to .D-dimensional de Sitter and carry out the 
field-dependent gauge transformation that enforces exact, Lorentz gauge. Of 
course this gauge transformation is ambiguous up to a homogeneous solution, 
which is uniquely determined by de Sitter invariance but which we leave un- 
specified at this stage. In section 4 we describe the unique, de Sitter invariant 
solution that was found by solving the propagator equation in Lorentz gauge 
[T9"] . In section 5 we show how the hitherto unspecified, homogeneous part 
of the gauge transformation from section 3 can be chosen to make the two 
propagators agree. Our conclusions comprise section 6. 

2 Deriving Average Gauges 

The purpose of this section is to explain how one derives average gauges 
from the exact gauges of the canonical formalism. We shall take flat space 
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quantum electrodynamics (QED) as the object of study and first sketch the 
technique for passing from Coulomb gauge to Feynman gauge on the manifold 
R D . We then consider the same process for the manifold T " 1 x R to show 
explicitly why full Feynman gauge cannot be imposed. 
The QED Lagrangian is, 



C = —F^F^ + ^(id^-eA^ - 



(25) 



The canonical dynamical variables of Coulomb gauge (V- A(t, x) = 0) are the 
transverse components of the vector potential, Af (t, x), and the electric field, 
Ef(t,x), as well as ip(t,x) and ip(t,x). The \i = component of the vector 
potential is not an independent dynamical variable but rather a nonlocal 
functional i[)](t,x) of the charged fields, 

1 



A = 



e^7 ijj 



(26) 



And the Hamiltonian density of Coulomb gauge is, 

H = \EjEj + -F^ - ~d> V 2 $ + + ej'Aj + m)V . (27) 

The usual connection between the canonical and functional formalisms is 
made through the vacuum expectation values of T*-ordered products of op- 
erators. (T*-ordering is the same as time-ordering except that derivatives are 
taken outside the ordering.) Suppose 0[ip, ip, E T , A T ] represents an arbitrary 
functional of the canonical operators. The fundamental relation between the 
canonical and functional integral formalisms is, 

T 



n 



n 



T*(o[ifj,tfj,E T ,A 7 

[chji] [dip] [dE T ] [dA T ] 



,i Jd D x[i^°i,+EfAf-H] Q 



V*j ^/S F T , A T 



.(28) 



Because they will play no role in our analysis we have suppressed the initial 
and final state wave functionals. 

We henceforth denote expression (|2"8"j) as (O). A 7-step process of func- 
tional manipulations carries it to Feynman gauge on R D : 

1. Perform the Gaussian integrals over the transverse electric field, 



[dE T ] jS^lEjAj-WEl] ^ £T A 



0'\^,^,A T }e i I dDx ^ A T A f , (29) 
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where (D'[ip, ip, A T ] is 0[ip, if), A T , A T ] plus the delta function correlator 
terms that arise from eliminating the various pairings of E T, s. 



2. Restore the longitudinal component of the vector potential, 



[dA 1 



[dA}8 



V-A 



det(-V 2 ) 



(30) 



Note that this gives the square root of the Faddeev-Popov determinant 
for Coulomb gauge. 



3. Restore the temporal component of the vector potential, 

e i/d^*V** / det( _ V 2) x 



[dA ] e { / d ° x ih^o^o-i'-r ^^} _ (31) 



Note that this gives the remaining bit of the Faddeev-Popov determi- 
nant. At this stage (O) takes the form, 



(O) 



[dip] [dip] [dA] 8 



V-A 



det(-V 2 



xe i fd D x[\A i A i +\d i A d i A -\F lj F ij +^{i^d tl -YeA^-m)^\ Q"'^^^ ? (32) 

where 0"[if), if), A] is 0'[if), if), A] with possible factors of $ [if), if)] replaced 
by A (ii desired; it is not necessary) and with the addition of appro- 
priate correlator terms for pairings of Ao's. 

4. Express the integrand as an invariant. Any good gauge can be used 
to express an arbitrary functional of the fields as a gauge invariant 
which happens to agree with the original functional when the gauge 
condition is obeyed [21 [20]. We do this for the action and for the 
operator 0"[ip,ip, A], 



V-A 



x e 
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i fd D x[\A i A i +\d i A a i A \ 



V-A 



x O" 



if),if),A 



V-A 
V-A 



x e 



fd D x±F oi F oi 



x a 



if), if;, A 



After invariantizing in this way our expression for (O) is, 



(O) 



[dip][dip}[dA}8 



V-A 



det(-V 2 ) e iS ™ O 



(33) 
(34) 



(35) 
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5. Make a functional change of variables to Lorentz gauge. Consider the 
field- dependent gauge transformation, 



4 = A,-d,9 1 [A}: (36) 
if/ = e ie9l[A] x ib , (37) 
i/ = e~ ieei[A] x ^ (38) 

where the gauge parameter is, 

Oi[A] = -^A . (39) 

Because S inv and O inv are gauge invariant, only the gauge fixing delta 
functional and the measure will change. To get them, note that the 
inverse transformation for the vector potential is, 

A, = ^ - d^A> . (40) 

It follows that the Coulomb gauge condition on implies the Lorentz 
gauge condition on A' , 

d,Ai = PA'^ . (41) 

And the functional Jacobian converts the Faddeev-Popov determinant 
to the one appropriate for Lorentz gauge, 

[dA] det(-V 2 ) = [dA'] det(-<9 2 ) . (42) 

We can therefore write (O) as, 



(O) 



[dif/] [dip'] [dA'} S [d^A'J det(-<9 2 ) e iS ™ O inv . (43) 



6. Add a inhomogeneous, (D -number term to the gauge fixing functional. 
Make an additional change of variable, 

K = (44) 

if>" = e ied2lf] x ifj , (45) 



—HI 



ib" = e - ie(>2[f] x ib , (46) 



9 



where the gauge parameter is defined in terms of an arbitrary C- number 
field /(*,£), 

«] = ^/- (47) 
After the transformation we have, 



(0) = 



[#"] [dip"] [dA"\ 5 [d^l - f] det (-<9 2 ) e lS '™ O mv . (48) 



7. Functionally average over the inhomogeneous term. By construction 
(O) has no dependence upon the function f(t,x). It is therefore un- 
changed if we multiply by a normalized Gaussian and functional inte- 
grate over /, 



(O) 



[df]e l I dDx ^ f2 x (O) , (49) 
"] [#"] [dA"\ det(-<9 2 ) e i5Fc " O im . (50) 



This is the Feynman gauge functional formalism we sought. 

Let us see how the sequence of functional manipulations described above 
changes when the noncompact spatial manifold is replaced by the com- 
pact manifold T D ~ l . A major difference is that Fourier integrals become dis- 
crete sums. Suppose the range of each spatial coordinate is — L < x l < +L 
for i = 1, 2, . . . , (D — 1). Then any function f(t,x) can be expressed as a 
discrete Fourier sum, 

f(t,x) = £ f n (t)e ika - s , (51) 

where k = n/L is the fundamental wave number. Note that the action of V 2 
on any such function annihilates the n = mode, 

V 2 /(t,f) = - £ (kn) 2 f n (t)e™ L - ln - s . (52) 

Hence the instantaneous Coulomb potential can only be defined for con- 
figurations of ip(t,x) and ip(t,x) which have zero total charge (this is the 
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linearization stability constraint!) and the resulting potential has no n = 
mode, 



${t,x) 



1 

V2 



(53) 



E 



e ikn-(x-s>) ^ f ') 7 o e ^,( t? ^) _ (54) 



Of course this means that when A (t, x) is restored in step 3, it cannot contain 
any n = mode. 

Another important change concerns restoring the longitudinal part of the 
vector potential in step 2. Under a gauge transformation the Fourier mode 



noninteracting theory they would behave like free quantum mechanical par- 
ticles rather than harmonic oscillators.) A second consequence is that the 
Coulomb gauge delta functional lacks a n = mode, 



To recapitulate, the changes associated with working on the compact 
spatial manifold T D_1 are: 

• The field A (t,x) contains no n = mode; and 

• The Coulomb gauge delta functional contains no n = mode. 

It is immediately obvious that, while we can still make the functional change 
of variables in step 7 to enforce exact Lorentz gauge, neither the resulting 
A' (t, x) nor the Lorentz gauge delta functional will contain an n = mode. 
It follows that there is no valid way to add the n = mode of the Feynman 
gauge fixing term. So the result is just what we expected: adding the full 
Feynman gauge fixing term is incorrect, and all conclusions drawn from this 
formalism are suspect. It should be clear that the same sort of problem must 
occur as well in de Sitter, and for gravitons as well as for electromagnetism, 
so we now have a proof that the average gauges for which de Sitter invariant 
solutions have been found are not valid. 





(56) 



tG-Rn^O 



11 



Let us return to the context of flat space electromagnetism on the man- 
ifold T D ~ X x R and consider what goes wrong if the problem we have just 
demonstrated is ignored and the covariant gauge fixing term (with n = 
mode) is erroneously added to the invariant action. In that case there is an 
extra, homogeneous contribution to Ao(t, x), which should not be present, 
and a corresponding extra contribution to propagator. For a point charge q 
this extra term produces a homogeneous contribution to A (t, x) which grows 



One might object that the extra term is harmless because it makes no con- 
tribution to the electric field, however, the undifferentiated potential does 
contribute to the interaction energy. A special case of some interest in quan- 
tum electrodynamics is the interaction of a particle with its own force fields. 
In this context it has been noted that using the de Sitter invariant, Feyn- 
man gauge propagator [21 J (which must also contain a spurious homogeneous 
part) results in on-shell singularities for the one loop self-mass-squared of a 
charged scalar [18]. Just as the analysis of this section suggests, these on-shell 
singularities disappear either: 

• When using a non-de Sitter invariant gauge on the open coordinate 
sub- manifold (which has no linearization instability) [IS] ; or 

• When using the de Sitter invariant Lorentz propagator (which is exact) 
[22]. 

It should be emphasized that there is no mistake in the Allen- Jacobson solu- 
tion for the Feynman gauge propagator [21] ; it is the gauge fixing functional 
which is at fault. 

3 Imposing Lorentz Gauge on de Sitter 

The purpose of this section is to make a field-dependent gauge transformation 
which carries the photon propagator from a non-de Sitter invariant average 
gauge, defined on the open coordinate submanifold, to the de Sitter invariant, 
exact Lorentz gauge, which can be extended to the entire de Sitter manifold. 
We begin by reviewing the geometry and coordinate system, then we give 
the noninvariant gauge condition and the associated propagator. The next 



like t 2 , 




(57) 
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step is making the transformation. Of course this is ambiguous up to surface 
terms which we leave to be specified in section 5. We close by decomposing 
the transformed propagator (without the homogeneous contributions) up into 
two convenient pieces. 

We work on the open conformal coordinate submanifold of D-dimensional 
de Sitter space. A spacetime point x^ can be decomposed into its temporal 
(x°) and spatial x % components which take values in the ranges, 

— oo < x° < and and — oo < x l < +00 . (58) 

In these coordinates the invariant element is, 

ds 2 = g tiv dx il dx v = a 2 rj PiU dx' l dx u , (59) 

where r]^ is the Lorentz metric and a x = —1/Hx° is the scale factor. The 
parameter H is known as the "Hubble constant" . 

Most of the various propagators between points x^ and can be ex- 
pressed in terms of the de Sitter length function y(x; z), 

2 /. n n. \2 



y{x; z) = 



x — z 



-(\x°-z°\-ie) . (60) 

Except for the factor of ie (whose purpose is to enforce Feynman bound- 
ary conditions) the function y(x; z) is closely related to the invariant length 
£(x; z) from x^ to 

y(x;z) =4 sin 2 (^H£(x;z)) . (61) 

Because y(x; z) is a de Sitter invariant, so too are covariant derivatives of it, 
dy(x; z) 



dx^ 



H 2 a x a z (y5l5» + 2a z H&xX-2aJlHAx u -2ri^) . (64) 



= Ha x (y5l + 2a z HAx,) , (62) 
= Ha z (y5° u -2a x HAx u ) , (63) 

d 2 y(x; z) 
dx^dz v 

Here and subsequently Ax^ = rj^ (x — z) v . 

Electromagnetism is conformally invariant in D = 4 dimensions, which 
means that it takes the same form in conformal coordinates as in flat space. 
This is obvious from the gauge invariant Lagrangian, 

Anv = -- A F, v F pa g»o 9 ™^ = -] a D -% u F paV ^ . (65) 
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The wonderful simplicity of using known flat space results will not be pre- 
served if one adds any multiple of the de Sitter invariant, Feynman gauge 
fixing functional, 



as 



i 



D-4 



fr u A^-(D-2)HaA 



(66) 



(A semi-colon denotes covariant differentiation whereas a comma stands for 
the ordinary derivative.) However, a very simple formalism results from 
replacing the factor of (D — 2) with (D— 4), 



£ Nd s = -\a D -\r,^d^A v - (D—A)HaA( 



(67) 



With this gauge fixing functional the propagator takes the form [TOl [18] , 



NdS 



(x; z) = a x a z iA B (x; z){r]^ + <5°<5°) - a x a z iA c (x; z)5°5 
where the de Sitter invariant scalar propagators are, 



0*0 

V 1 



(68) 



iA B (x; z) 



iA c (x; z) 



B(y(x; 
C[y(x; z)j 



H ~ D T{D nS MD-2,i;-;i-y 

(4tt)t T(f) H 



D 
~2 



(4tt)t r 



2 / 



(69) 



HD '^ D -\ Fl (D-^-y) . (70) 



One nice thing about (168]) is that its tensor factors are constants. Another 
is that each of the scalar propagators which multiply them consists of the 
conformal propagator plus a series of less singular terms which vanish in 
D = 4 dimensions, 



B(y) 



H 



D-2 



(47f)^ 



C(y) 



H 



D-2 



(47T 



D 



v 2 



4xf-i 



4xf-i 



E 

n=0 L 



E 

n=0 



r(n+f),y 



r(n + 2) M 
T(n + D-2) ,y 



r(n+f; 



(71) 



r(n+f-l) m n-f+2 



(I)" 



(n+1) 



T(n+2) V4 

r(n-KD-3) /j/N 
T(n+#) W 



(72) 
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So the infinite series terms only need to be retained when they multiply a 
potentially divergent quantity. Because the higher values of n vanish more 
and more rapidly at coincidence (that is, for y — 0) only a finite number of 
these extra terms ever needs to be included. 

It is now time to make the field-dependent transformation to Lorentz 
gauge, 

A'^x) = A,(x) - d,6[A)(x) . (73) 

This would be step 5 in the scheme of the previous section. The fact that 
A'^ obeys Lorentz gauge implies a differential equation for 9[A], 



(74) 



Of course there are many solutions, related to one another by homogeneous 
terms. Any choice of homogeneous term will enforce Lorentz gauge, whereas 
there can be at most one choice which gives a de Sitter invariant propagator. 
Because Section 5 is devoted to establishing de Sitter invariance and cor- 
respondence with the known solution [19] , we postpone specification of the 
homogeneous term until then. For now we express the solution in a general 
way, 

' d D x' G{x; x') 



6[A\{x) 



d 



-g(x')g^(x')A a (x') 



(75) 



iv ' dx' p 

where G(x; x') is some Green's function of the scalar d'Alembertian, which 
we specify in the next section, and V is some region of the manifold. The 
actual solution for 6*L4](x) consists of (1751) — with definite choices for G(x; x') 
and V — plus a definite homogeneous solution. For now we study the field 
transformed with only #[A](a;), 



A^x) = A^x) 



dJ[A](x) 



(76) 



The transformed propagator is the vacuum expectation value of the T*- 
ordered product of two A's. Because T*-ordering moves any derivatives out- 
side the time-ordering symbol we can express this as, 



n 







dx^ Jv 




A^x)A u (z) 

d D x' G(x; x) x 

d D z'G{z;z') x 



(p 


T 


d 




dx'P 


d 




dz' a 





A^{x)A v {z) 



g(x')g p °(x')(nT A a {x')A v {z) 



-g(z')g a ^z')(QT A,(x)A p (z> 
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X 



d d 
dx'p dz' a 



-g(x')g^(x')J-g(z')g^(z')(Q T A a (x')A p (z>) 



• (77) 



By substituting the noninvariant propagator (1681) and using the fact that 
y(x; z) depends upon spatial coordinates only through their difference, we 
can write the three differentiated, square-bracketed terms as, 



d 

dx'p 



_d_ 

dz 71 



-g(x')g^(x'){nT[A a (x')A u (z) 
+t°»{^[^- 1 a z B(y(x'-,z) 



-g(z')g^(z')(QT[A,(x)A^z r 
d 



+5 "\ Or" 



Oi X Oi z 



n 



+ 



n 



+ 



dz 



d 



dx'° 



d 
dz 1 ® 



d 8 
dx'p dz' a 
d d 



dx n dz' 1 



-g{x')g^{xy-g{z')g a %z')(nT[AM)Mz') 

d d 



d 

Y a®r l a z B(y(x'] z)^ 

a^r 1 a z C(y(x';z))}^ (78) 
d 

--Q^[ a * a ° rlB (y^ z '))_ 

'a x afr 1 C(y(ar,z , ))]\, (79) 



,a T ia z > 



B[y{x ; z 



dx'° dz'° 



{a xl a z ,) D - l C{y{x'-z'))\ . (80) 



All of this suggests that we would do well to organize the transformed propa- 
gator into a doubly differentiated "Integral Term" and the remaining, "Other 
Term" , 



n 



A^x)A v (z) 



n 



-T(x; z) + Lo v ] (x; z) . (81) 



dxP- dz v 



The Integral Term is, 

I(x; z) = 

a z f d D x G^x^a^B^x'-z)) + a x f d D z'G(z; z^afr 1 B(y(x; z'fj 



v 

d d 



dx'° dz'° 



(a x >a z/ ) D - 1 B(y(x';z') 



{a^a^C^yix'^'))] . (82) 
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And the Other Term contains everything else, 



(x; z) = a x a z B(y(x; zfj [?7^+5j5°] - a x a z C(y(x; zfjd^l 



-5 



d 



u dx^ Jv 
d 



d D x'G(x;x')< 



d 



dz° 



dx'° 



a°r l a z B(y(x>-z)) 
a^ 1 a z C(y(x';z))}\-5^ J y d D z' G(z; z') 



x ■ 



'_d_ 

dx° 



a x afr 1 B(y(x;z'))}+ 7 ^[a x af,- 1 C(y(x; z'))]\. (83) 



dz'° 



It remains to act the derivatives to simplify our expressions for T(x; z) 
and [^O v ]{x] z). This is facilitated by some important identities obeyed by 
any function of y(x; z), 



n x F(y) 

d 2 F(y) 
dx°dz° 

d 2 F(y) 

dx l dz l 



iAir? 5 D {x — z) 
r(f-l) H D ~ 2 a D 

iArr^5 D (x — z) 
-1) {Ha) D - 2 



x Res[F] + H 2 UAy-y 2 )F" + D(2-y)F'\, (84) 



■D 



x Res[F] + a x a z H 2 < \8 — (Ay—y 2 )]F" 



H 



d 



dz° 



-{2-y)F' + p + ^i) [-2(2-y)F"+2F' 

a z a x 

a x a z H 2 L(2-y)F" - 2{D-l)F> - a(^ + ^)f") , 
^^]^)=^ 2 {-2(2-^ + 2(J + J)F'}. 



(85) 
(86) 
(87) 



Here □ is the covariant scalar d'Alembertian and Res[F] is the coefficient of 
in the Laurent expansion of F(y). We shall also require some identities 
specific to B(y) and C(y), 



(Ay-y 2 )B"(y) + D(2-y)B'(y) = (D-2)B(y) , 
(Ay-y 2 )C"(y) + D(2-y)C'(y) = 2(D-3)C(y) . 

And there is a very useful relation between B{y) and C(y), 



C(y) = \(2-y)B(y) + J^ 



where k 



(4vr)f T(f) 



(88) 
(89) 



(90) 
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Note finally that substituting (f9~0j) into (1591) and using (18"8|) implies, 

(^-y 2 )^) + (D-2)(2-y)B(y) = -2k . (91) 

It is best to start with the Other Term because it involves only first 
derivatives. The reduction is straightforward for the B term, 



_d_ 



a°- l a z B(y( 



x: z 



a? 1 a, 



a? a z H{ 2B' + — 



Ha zB + $L b' 



dz° 



I "2 



[2-y)B' + B 



(92) 
(93) 



We begin the same way with the C term but then use (1901 to convert most 
of the C"s to .B's and simplify with (}9Tj) . 



d 
dx° 



af 1 a 2 C(y(x; 2;) J 



a2 V 



(D-l)Ha x C + 



dx° 



a»a z H{ -{2-y)C + (D-1)C + — 2C' 



(94) 
(95) 



a°a z H{ 2C - ~{2-yfB' + UD-2){2-y)B + k + 



1 



1 



a z r 
a 



a z v, 



a°a z H{ 2C - 2B' + — {2-y)B'-B 



(2-y)B'-fi]| , (96) 
(97) 



Hence (I9"3"|) and (19~T|) almost completely cancel and our final result for the 
Other Term is, 



[x; z) = a x a z B(y(x; zfjrj^ + a x a z f^~y(x; z)B(y(x; zfj 
d 



k 



D-3 



-2HaJ 



d D x'^J-g(x')G(x; x')C(y(x'; z) 
, d 



d D z'^g(z')G(z;z')C(y(x;z')). (98) 



dz u Jv 

The Integral term T(x; z) involves second derivatives. We only need (1561 
to reduce the spatial case, 



d d 



dx l dz l 



(a x a z ) D l B[y(x]z) 
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'flu a. 



H\a x a z ) D \ 4(2-y)B" - 2(D-1)B> + U + 



-45" 



>. (99) 



Reducing the temporal derivative term is more involved. We begin by passing 
the derivatives through the scale factors, then employ relations (185]) and 
f lH7|) . and convert most of the C(y)'s to B(y) using fl90l) . eliminating second 
derivatives with (I88H89P as needed. The result is, 



d d 



dx° dz° 



{a x a z ) D - l C[y{x-z))] = (a^/" 1 

5 



f g g 

\ dx° dz° 



+(D-1)H 







+a. 



(D-lYa x a z H 2 }C(y(x;z)) ,(100) 



dz° dx° - 

ia°5 D (x-z') + H 2 (a x a z ) D \4C + 4(2-y)B" - 2(D + 3)B' 



• a z a x 

Adding the two terms gives a compact form 



-4B" + 2(2-y)B' -2B I .(101) 



d_d_ 

dx l dz l 



^cl x cl z j 



\D-1 



B(y) 



g g 
<9x° (92° 



\-D-l 



C(y) 



ia®5 D (x — z) 



+ (a x a z ) D H 2 l-AC + 8B' 



a, a x 



-2{2-y)B' + 2B [.(102) 



A further simplification can be effected by means of the identity, 



a x a z ) D n x —B(y) = ia»5 D {x-z) + {a x a z ) u H 2 { -AB'+^ 2{2-y)B'-2B 



,D X D 



Dtt2. 



(103) 



Using this and the result with x^ and z^ interchanged gives, 



d d 

dx l dz % 



\a x a. 



B{y) 



d d 
dx° dz° 



(a x a z ) D - l C(y) =ia u J u {x-z) 



D X D< 



+ (a x a z ) D \-AH 2 C-n 



—B 

a* 



a x 



>. (104) 



Our final result for the Integral Term involves the surface integral, 



Sr(x: z) 



-B(y(x; z) \ — j d D x \J —g(x') G(x; x r )n x < 
z JV 



a,./ 



B[y(x'; z] 



dV ^ v 



— iA B (x f ; z)d' u G{x- x')-G{x; x')d' v 

a 2 



— iAnix': z, 



,(105) 
.(106) 
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This function is obviously homogeneous; that is, \3 X annihilates it. In the 
next section we will show how to choose the homogeneous contribution to 
the full gauge parameter #L4](:r) so as to cancel it and similar terms. The 
final result for the Integral Term is, 



X[x\z) = I d D x' J —g(x') G(x; x')Sb{z; x') 

J V" 

+ J^d D z' ^ -g(z') G{z; z')S B {x; z') + i j d D x' j -g{x') G(x; x')G{z; x') 
-AH 2 j d D x'^/-g{x')G{x;x') j d D z'^-g(z') G(z; z')C(y{x'; z')) .(107) 



4 The Invariant Propagator 

The purpose of this section is to show that the transformed propagator of 
the previous section agrees, up to surface terms, with the unique de Sitter 
invariant propagator which was found by solving the propagator equation in 
Lorentz gauge [T9] . Of course we begin by describing that solution. We then 
decompose it in analogy with the scheme (|81[) of the previous section, into 
an "Integral Term" and an "Other Term." At this stage there is a digression 
to derive an identity for the convolution of scalar propagators. The section 
closes by applying this identity to demonstrate that the two propagators 
agree up to surface integrals. 

The Lorentz gauge propagator equation has a unique de Sitter invariant 
solution which can expressed in terms of a function 7(2/) [19], 



as 



1 



(x: z 



j d 2 y(x;z) 

A(D-1)H 2 \ dx»dz v 



-(%-y 2 ) 7 '-p-l)(2-y)7 
dy dy 



dx^ dz v 



( 2 - y ) 7 '-( J D-l) 7 .(108) 



The function 7(7/) has a very complicated series expansion, 

(n+l)T(n+D-l) 



1 /D — l\ H D ~ 2 
7(y) = o(nZa) 



2 \D-3' (4 



7T 



;{( C -3)r(f-i)(i)-' + f;" 

2 I Z if n=0 . 



X 



r(n+f +1) 

^(2-^)-^(^-l)+^(n + J D-l)-^(n + 2)](|) n 



(n-f + 3)r(n + f + 1) 
T(n+3) 



2 

D 
2" 



D 
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+^(™ + y + l)-^-y+4)](f) 



4 



(109) 



Although it might seem unwieldy, this formalism has been used to perform 
several two loop computations in scalar quantum electrodynamics [22| [23]. 
The function 7(2/) obeys the second order differential equation, 



(42/-y 2 ) T " + {D + 2){2-y)i - 2(D-1) 7 = 2(D-l)B'(y) 
One consequence is, 

d ' 



;no) 



<9y 



( 4y _^)y_( jD _l)(2_ 2/ ) 7+ 2( J D-l)S 



{2-y)i-{D-l) 1 . (Ill) 



Hence we can decompose the invariant propagator in a form analogous to 
that of the transformed propagator (IHTT) . 

d 2 y(x; z) 



A dS 



[X] Z) 
1 



2# 2 

9 9 



B[y(x;z) 
I 



dx^dz v 

-(Ay-y 2 ) 1 '-(D-l)(2-y) 1 + 2(D-l)B 



■(112) 



4(D-l)H 2 dx» dz v 
where the notation "/[/]" of a function /(y) stands for its indefinite integral, 

I\f\(y) = fdy'fW)- (H3) 

At this point it is useful to digress on the subject of scalar propagators. 

Three were introduced in the previous section H x G(x; z), iAb{x\ z) and 

iAc(x; z) — and it might seem that there is a bewildering variety of them, 
each with its own important special properties. However, a unified treatment 
can be given in terms of the equation, 



-g(x) I D x + 



D - 1n2 



H 2 UA„(z; z) = i5 D (x-z) . 



(114) 



The three propagators of the previous section correspond to the following 
choices for v: 



i x G(x; z) 
iA B (x; z) 
iA c (x; z) 



v 



v 



D — 


1 


2 




D — 


3 


2 




D — 


5 



(115) 
(116) 
(117) 
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For general v the spatial plane wave mode functions corresponding to Bunch- 
Davies vacuum are, 



u v (x°,k) = a-^H?\-kx°). 

When it exists, the Fourier mode sum for the propagator is (2 
d D - x k 



iA u (x; z) 



(2n] 



D-l 



118) 



e{x {3 -z [) )u v {x\k)u* u {z^k) 

+6(z -x )u* u (x , k)u u (z°, k) 1.(119) 



When this sum exists the result is de Sitter invariant [25 

fjD-2 p/- -D-l i .AW-D-l 



iAJx\ z) 



(47t)t 



r(f; 



„ ,D- 1 D-l D v 
FA +v. v. — ; l-~ 

l \ 9, '2 ' 2 ' 



4' 
(120) 

When the Fourier mode sum f)119p is infrared divergent it must be cut 
off, either by making the mode functions less singular for super-horizon wave 
lengths [26J or by working on a spatially compact manifold [27]. Either 
procedure breaks de Sitter invariance. A special case of some importance to 
our discussion is v = (D — l)/2, for which the result is [28l 124] . 



iA A (x; z) = A(y(x; zn + k\n( 



;i2i) 



where the constant k was defined in ( 1901) and the function A(y) has the 
expansion, 



H 



D-2 



(4vr)i 

oo 

+ E 

71=1 



1 r(w+D-l) /y 
n T(n+#) U 



4 N 2-1 r(# + l) /4n#-2 



t-2 v 2/ 



— 7T COt — — 



TTDxTp-l) 



v 2 / r(f 



a) 



1 r(n+f+l)/y 



n-f + 2 T(n + 2) 



■(f) 



.(122) 



As with the expansions (I71H72P for B(y) and C(y), the infinite series terms 
of A(y) vanish for D = 4, so they only need to be retained when multiplying 
a potentially divergent quantity, and even then one only needs to include a 
handful of them. This makes loop computations manageable. For a massless, 
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minimally coupled scalar with a quartic self-interaction, two loop results have 
been obtained for the expectation value of the stress tensor [25] , for the scalar 
self-mass-squared [29] and for the quantum-corrected mode functions [30J. 
In Yukawa theory it has been used to compute the expectation vlaue of the 
coincident vertex function at two loop order [31], and it has been used for a 
variety of two loop computations in scalar quantum electrodynamics [22, 23J. 
It should also be noted that the de Sitter breaking correction to zA^(a;;z) 
in expression (11211) can be derived from the infrared-truncated mode sum 
|24j . and it serves to reproduce the classic and well known result for the 
coincidence limit of the propagator |32j . 

The function A(y) obeys a differential equation analogous to (1551) and 



(4y-y 2 )A" + D(2-y)A' = (D-l)k . (123) 
A number of identities relate the derivative of A(y) to B(y) and C(y), 

1 



A' 

(4y-y 2 )A' 



■(D-3)B + C, 
-2(D-2)B - (2-y)k . 



(124) 
(125) 



It is also useful to note the result of acting the scalar d'Alembertian on a 
function of the scale factor, 



□/(a) 



a 2 f"(a) + Daf'(a) 



(126) 



Now consider Green's second identity for any two functions F(x') and 
G(x>), 



F{x')^-g{x') n'G(x') - G(x')^/-g(x') U'F{x') 

~F{x')d' v G{x') - G(x')dlF(x')} 1.(127) 



= d'^-g(x')g^(x' ^ 

We choose G(x') to be any symmetric Green's function G(x; x') = G(x'; x 

-g(x) UG(x;x') = 5 D (x-x' 



(128) 

We can obviously integrate (I127P over any region V with boundary dV to 
conclude, 



F(x) = fd D x'^-g(x') G{x; x')\3 x ,F(x') 

F{x')d' u G{x; x') - G{x; x')& v F{x') 



d D -y, 



-g'g^ 



■ (129) 
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Relation (I129p is true for any Green's function so we are free to use the A- 
type propagator, G(x;x') = — i x iA A (x;x'). Relation (11291) is also valid for 
any function F(x) so we are free to make the choice, 



F(x) 



iA u (x; z) —iA A (x; z) 

[ (i2_l)2_ z/ 2] H 2 



(130) 



The surface terms involving iA A obviously cancel so the result is, 

iA u (x; z) — iA A (x; z) 

[(^_l)2_ z/ 2] i/ 2 



i d D x'J-g(x')iA A (x;x')iA 1/ (x';z) 



+i d D -y^-g'g^ 

JdV H v 



iA v (x'; z)d' a iA A (x; x')-iA A (x; x')d' a iA v (x'] z) 

[(^ r i)2- Z y2] i J2 



.(131) 



We call fllUTj) the "Convolution Identity" 

Choosing v = (D — 5)/2 in the Convolution Identity f 1 1 3 1 j) gives us a 
relation for the C-type propagator C(y), 

iAcix: z)-iA A (x: z)-Sn(x: z) 



i d x y — g(x') iA A (x; x')iA c (x'; z) 



2(D-3)H 2 



(132) 



where the surface term is, 



S c (x;z) = d D -yJ-g>g'<" iA c {x'- z)& a G{x; x')-G{x- x')& a iA c {x'- z) . 

JdV ^ v L J 

(133) 

We now substitute (I132p in our result fl98l) for the the Other Term of the 
previous section, 



(x; z) = a x a 



■B(y(x; z))^ + a x a z S^y(x; z)B( 



y{x;z) 



k 



-2H 



aJl dS c (x;z) 



D-3, 

iA c (x; z)-iA A (x] z) 
2(D-3)H 2 
aJl dS c (z;x) 



■ 5°S° 

w fJ, w V 



+ 



(D-3)H dx» (D-3)H dz> 
The derivative is easy to simplify using (11241) . 



.(134) 



o d \ iA c (x;z)-iA A (x;z) 



2(D-3)H 2 
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(D-3)H dx» 



a - a') + 



— 5y u a x a z , (135) 



= a x a z {-a z HAxXB - \yB&$, + <^J.(136) 

Combining everything results in an expression for the Other Term which is 
almost de Sitter invariant, modulo the surface terms, 



(x; z) = a x a z < --ybffi - a z HAx^ v + a x 5^HAx v + rj^ >B 



D-3 



a x a z b„b u + 



H v 



a z 5° dS c (x; z) 
(D-3)H dx* 



+ 



aJl DSo(«x) 



(D-3)H dz" 



d 2 y(x;z) B ka x a z 5ffi a z S° u dS c (x;z) a x 6° dS c (z;x) 
dx»dz v 2H 2 D-3 D-3 dHx» D-3 dHz u ' [ ' 

The first term on the right hand side of (11381) agrees with the first term 
in our decomposition (11121) for the invariant propagator. We must obviously 
choose the homogeneous contributions to the gauge parameter 6 l [y4](x) so as 
to cancel the surface terms in (I138p . That leaves the term proportional to k, 
which relation (11261) allows us to recognize as a potential part of the Integral 
Term, 



ka x a z 5®5® 



d d [ k In (a r a. 



+ const x ln(a x a z 



(139) 



D-3 dx» dz v \2(D-3)H 2 } 

We will presently see that precisely the bracketed expression is needed to 
make 1{x\ z) de Sitter invariant, up to surface terms. 

A matter of great importance for us is what the Convolution Identity 
(11311) gives when the index v is chosen to be (D — 1)/2, corresponding to 
the A-type propagator. The term on the right hand side obviously gives a 
derivative with respect to the index v, 



lim 



iA u (x; z) —iAa{x; z) 

[(^_i)2_ z/ 2] i/ 2 



§^iA u (x; z) 



(D-1)H 2 

_ iAa'(x;z) 
~ {D-l)H 2 ' 

Hence the convolution of two A-type propagators gives, 

iAA'ix; z) 



i d D x' \J —g{x') iAa{x] x')iAa{x'] z) 



■ + 



S A {x\ z) 



(D-1)H 2 (D-l)H 2 



(140) 



(141) 



(142) 
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where the surface term is 
S A (x; z) = ld 



p 

dV 



g'g' pa 



iA A ,(x'; z)d f a G(x; x')-G(x; x')&JA a ,(x'] z) 



(143) 

Like the A propagator, the A 1 propagator breaks de Sitter invariance. 
The simplest way to see this is by differentiating relation fjl 14|) with respect 
to v and then setting v — (D — l)/2, 



dv 



V 



D-U2 



H 2 



iAJx; z) 



U x iA A ,(x; z) + (D-l)H 2 iA A (x;z) . 



(144) 
(145) 



Because iA A (x; z) = A(y) + k\n(a x a z ) has a de Sitter breaking term it is 
clear that iA A > (x; z) must as well. From relation f)126p we infer, 



iA A '(x;z) = A(y(x\ z)) + const x \ia(a x a z ) + -k\ia 2 (a x c 



(146) 



We do not require the coefficient of the ln(a x a z ) term but the series expansion 
for A(y) is, 



Av) 



H 



D-2 



(4tt)t 



4x^-2 



oo 

(47t)t n=0 



n-f + 3 



r(w + f + 2) 
(n+2)\ 



^(2-^)-^(f-l)+^(n+y+2)-V( 



D 

n-- + 3 



\y) n+1 T(n + D) 



n+l T(n+f + l) 



^(2-^)-^(^-l)+^(n+D)-^(n+l) 



,(147) 



One can hardly fail to notice the similarities in the series expansion (11091) 
for j(y) and the expansion (I147P we have just given for the de Sitter invariant 
part of the A' propagator. The relation between them is, 

AV) = \(^)l[^y-y 2 W + (D-l)(2-y) 1 }-^(D-2)I[B]+const . (148) 

It is tedious but straightforward to check f!148j) using the series expansions 
but a simpler way of recognizing it is to act the scalar d'Alembertian on both 
sides. In view of (11451) the left hand side gives, 



j^A(y) = -(D-l)A(y) + const . 



(149) 
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To compute the right hand side we need the lemma, 

(2-y)(Ay-y 2 )i + (Ay~y 2 ) 1 + D(2-y) 2 1 = 2(D-l)l[(2-y)B'\ +const . (150) 

This follows from differentiation with respect to y and using the equation 
( 111 Oft for j(y). Now act D/H 2 on the first term on the right hand side of 
(11481) . then use the 7 equation (IllOp . and finally relations (19~T1) and (11501) . 



(%-y 2 ) 7 '+0D-i)(2-y)7] = (4y-y 2 )[(4y-y 2 ) 7 // +p + l)(2- 2 /) 7 / 
-p-l) T ] + D (2-y)[(Ay-y 2 )i + (D-l)(2-y) 1 ] ,(151) 

(D-l)\(2-y)(4y-y 2 )i + (4y-y 2 ) 1 + D(2-y) 2 1 + 2(4y-y 2 )B').(l52) 



2{D-1)\(D-1)I {2-y)B' -(D-2){2-y)B'+const\. 



(153) 



Acting on the right hand side of (I148P and using identities (1911 . (190]) and 
( 11241) proves the relation, 

^{li^^^y-y'h'HD-m-yh} - 1 -{d-2)i[b\\ 

= ~(D-3)(D-l)l[(2-y)B'] - ~(D-3)(D-2)(2-y)B 

-(D-2)(2-y)B + const , (154) 

(155) 
(156) 



l(D-3)(D-l)I[B] - l -{D-l)(2-y)B + const 
-(D-1)A + const . 



The point of enduring all this analysis is that we can now recognize the 
"Integral Term" of the invariant propagator (11121) as a collection of propa- 
gators plus the In 2 (a x a z ) term of expression (I139p . 



X(x; z) 



A(D-1)H 2 
A 

(D-3)H 2 
iAa.'(x] z) 
~(D-3)H 2 



-{Ay-y 2 )i -{D-l){2-y) 1 + 2{D-l)B , (157) 
C—A 

(L>-3) 2 ff 2 +const ' ( 158 ) 
iA c (x; z)-iA A (x;z) 



(D-3) 2 H 2 



+ 



A; In (a x a z ) 
2(D-3)H 2 



+const x \ia(a x a z ) + const .(159) 
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Note that the two unknown constants are irrelevant because they drop out 
when one differentiates with respect to x M and z v to get the propagator. 

We can make contact between the Integral Term (11591) of the invariant 
propagator and the Integral Term fjlOTj) of the transformed propagator by 
expressing the propagators as convolution integrals, 



iA A /(x; z) 



D-l 



d D x'J-g(x') G(x; x')iA A (x'; z) 



(D-3)H 2 \D-VJv 
iA c (x; z)-iA A (x; z) 
(D-3) 2 H 2 
S c (x; z) 



S A {x]z) 
(D-3)H 2 



,(160) 



d D z'J-g{z')G{z-z') 



~i. d x \/-9(x')G(x; x') 
iA c (x;z') 



, %A c {x'\ z) 



D-3 



S c (z;x) 



2{D-3) 2 H 2 Jv 
Now break up the prefactor of (I160p as, 

sD-U 1 



D-3 



1 



D-3 2(D-3) 2 H 2 
1 



.(161) 



(162) 



D-3 D-3 

and combine the convolutions multiplying the last two factors with the convo- 
lutions of (11611) to produce the combination iAq — iA A that can be recognized 
as another convolution, 



iA A /(x; z) 
(D-3)H 2 



iA c (x; z)-iA A (x; z) 



(D-3) 2 H 2 
ij^ d D x 1 \J —g{x') iA A (x; x')iA A (x'; z) 

iA c (z;x')-iA A (z;x') 



v 



+i / d D x'J-g(x') iA A (x; x') 



+iJ v d D z'^-g(z>)iA A (z ] z') 
S A (x; z) 



D-3 

iA c (x; z')-iA A (x; z') 



D-3 
S c (x;z) 



S c (z;x) 



(D-3)H 2 2(D-3) 2 H 2 2(D-3) 2 H 2 



,(163) 



-i d D x \J — g(x') iA A (x; x')iA A (x'; z) 

+AH 2 J v d D x'^-g(x') iA A (x; x')Jd D z yj-g(z') tA A (z; z') iA c (x'; z') 



r v d D x'^-g(x')G(x;x')'- 



<Sc(z;x') 
D-3 



d D z^-g(z')G(z;z')- 



S c (x\ z') 



D-3 



Sa(x; z) 



S c {x;z) 



S c (z;x) 



(D-3)H 2 2(D-3) 2 H 2 2(D-3) 2 H 2 
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.(164) 



We obtain the desired relation by adding the Integral Term (I107p of the 
transformed propagator to the In 2 (a x a z ) contribution 01391) from the Other 
Term (and some pieces which drop out when differentiated by x M and z u ), 



Z(x; z) + 



k In 2 (a x a z ) 
2(D-3)H 2 

= T(x; z) 



+ const x \n(a x a z ) + const 



Sa(x; z) 
(D-3)H 2 



Sc(x:z)+S c (z:x 



2(D-3) 2 H 2 



+[ d D x'^J-g(x')G(x;x') 



Sn(z\ x') 



+ d D z'^-g(z')G(z;z') 



S B (x: z)-\ 



<Scjz;x') 
D-3 



^ i S c (x; z') 
D-3 



.(165) 



5 Determining the Homogeneous Part 

The purpose of this section is to show that we can make the transformed 
propagator agree with the invariant one by correctly choosing the homoge- 
neous part of the full gauge parameter #L4](:r). We begin by summarizing 
the relevant results of the previous two sections concerning the gauge param- 
eter given in ( 17511 . which enforces Lorentz gauge but not de Sitter 
invariance. The resulting propagator agrees with the invariant one ( 1108111091 ) 
up to three surface terms which we denote as M-type," "5-type" and "C- 
type" according to the mode functions which they involve. We then exhibit 
a homogeneous gauge parameter A0L4](x), depending upon Aq, which can 
be used to absorb the B-tjpe and C-type surface terms. The section closes 
by deriving a homogeneous gauge parameter 56 l [A](a;), depending upon Ai, 
which absorbs the A-type surface terms and results in complete agreement 
with the invariant propagator. 



5.1 Summary of Previous Results 

Our goal is to construct a functional change of variables that is also a gauge 
transformation, 

M l Sx)=A li {x)-d lt 9[A}{x). (166) 
We want the field-dependent gauge parameter 0LA](x) to do two things: 

1. Make the field A' (x) obey Lorentz gauge; and 
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2. Make the propagator associated with A' (x) agree with the unique, de 
Sitter invariant solution of the Lorentz gauge propagator equation 

The first condition implies a second order differential equation for 0[A](x) 



-gU9 = dp y/=ggTA 



(167) 



Of course this only defines 0L4](x) up to a term which is annihilated by 
the scalar d'Alembertian. Because propagators obey Feynman boundary 
conditions we took the inhomogeneous solution to be the convolution of — i 
times the scalar propagator (G(x; z) = —% x iAA_(x; z)) with the right hand 
side of ([T67D , 



9[A]{x) 



d D x' G(x; x') 







Ox'p 



-g(x')g^(x')A a (x') 



(168) 



The result of just performing this transformation defines a field, 

A(x) = A^x) - d$[A\{x) , (169) 

which obeys the Lorentz gauge condition but whose propagator does not 
quite agree with the invariant one. 

We decomposed the propagator of A^{x) into the double gradient of an 
"Integral Term" ffTOTI) and an "Other Term" ([98]), 



A^(x)A v (z) 



n 



d d 



l(x; z) + 



(x; z) 



dx^ dz v 

The invariant propagator can be broken up in similar fashion ( 1 1 12[) . 



(170) 



as 



d d 



[x: z 



l(x: z) 



B[y(x;z] 



d 2 y(x; z) 



dx»dz v y ' 2H 2 \ y dx»dz v 
It is desirable to shift the double gradient of a spatially constant term 

k In 2 (a T a? 



AJ(x; z) 



2(D-3) 



const x \n(a x a z ) + const 



(171) 



(172) 



from [ At C jy ](x; z) to 2(x; z). When this is done, the difference between what 
we want the full transformation to produce and what the 6*[y4](x) transfor- 
mation actually gives is, 



d 2 Al 

dx^dz v 



Z S° dS c (x; z) aj° dS c (z;x 



D-3 dHx» 



D-3 dHz" 



(173) 
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I(x; z) — 1{x\ z) — AX(x; z) 



<Sa(x;z) 
(D-3)H 2 



S c (x; z)+S c (z;x 
2(D-3) 2 H 2 



-I d D x'J-g(x')G(x;x') 



S B {z; x' 



S c (z;x') 



- d D z'^-g(z')G(z;z') 



D-3 
S B (x; z')-{ 



S c (x;z') 



D-3 



.(174) 



Each of the surface integrals, Sf{x; z), consists of a Dirichlet and a Neumann 
contribution, 



s F (x-z) = id D -y 

av f 



-g'g' pa 



F(x';z)d' a G(x;x')-G(x;x')d' (7 F(x';z) . (175) 



The functions F(x; z) associated with the three integrals are, 



Sa {x; z) =>> F(x'; z) = iAA'(x'; z) , (176) 

Sb(x; z) =>> F{x';z) = — iAb{x']z) , (177) 

a z 

S c (x;z) =s> F(x'; z) = iA c (x';z) . (178) 

Note that each Sf{x;z) is homogeneous on the first argument x M . The in- 
tegral Sa{x'i z) is homogeneous on z^ as well, and also symmetric under 
interchange of x^ and z^ 1 . 



5.2 Absorbing the 5-type and C-type Surface Terms 



Rather than absorb all the surface terms at once it is simpler to first cancel 
those of the "Other Term," which must also reduce those that remain in 
the "Integral Term" to pure A-type (homogeneous on both x^ and z^\ We 
accordingly seek a homogeneous gauge parameter A#[A](a;) which cancels 
(11731) . Because this will also change the "Integral Term" we write out the 
full transformed field, 



M x ) 



A^x) 



d 
dx^ 



A9[A}(x) 



(179) 



The propagator of A is, 



n 



A„(x)A u (z) 
d 



n 



dx^ 



n 



d 



A6(x)A u (z 
d 



T 

Q) - 



A„(x)A u (z) 
d 



n 



dz h 



n 



A,(x)A9(z) 



n 



n 



A9(x)9(z)+9(x)A9(z)+A9(x)A9(z) Q) .(180) 
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The terms on the final line of (11801) must belong to the Integral Term ( I174p . 
and most of the middle line of (11801) must similarly belong to the Other Term 
(11731) . If we assume A0L4](a;) depends only upon A then the break is clean 
and we have, 



n 



a z S c (x; z) 



A0( lMoWj . v {D _ m 

A6(x)6(z) + 6(x)A6(z) + A6(x)A6(z)}\tt) = (l-type Terms 

S c {x]z)+S c {z]x 
2(D-3) 2 H 2 



;i81) 



-/ d D x'J-g{x')G{x-x') 



v 



S B (z: x') + 



S c (z;x') 



D-3 



r v d D z'^-g(z>)G(z;z>) 



S B (x: z') + 



S c (x;z') 



D-3 



.(182) 



It is straightforward to see that relation (11811) fixes the homogeneous part 
of the gauge parameter to be, 



-g'g"" 



Aq(x') , . , , Aq(x') 

o a G(x; x )— G(x; x )d a 



Combining (11681) and (11831) gives, 



(183) 



n 



A6{x)6{z) 



n 



(D-3)HJv 



d D z'G{z-z') 



m 9 



dz'° 



aB 1 Sc(x:z') 



(184) 

The surface integral Sc{x; z') has the form (I175p with the function F(x'; z') = 
iA c (x'; z'). Multiplying this by the factor of a®,~ 1 and taking the derivative 
gives an expression which we can simplify using relations (11241) and fl90l) . 



d 
dz^ 



afrHA c (x'; z') = Ha°{ {D-1)C - {2-y)C + 2^- C 



a,' 



(185) 



Ha" I 2C + (D-3)C - -(D-3)(2-y)B 



+ (£>_ 3 )^ B - (2-y)A' + 2— A' >, (186) 



Ha?,{2C+(D-3)^B\+a° 



-i d 



a z > 



dz'° 



iA A (x':z') 



187) 
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The final term involving iA^(x'; z') gives rise to an A-type surface term whose 
form we will work out in the next subsection. We can therefore write, 



n 



AB{x)B{z) 



n 



A— type Terms 



d D z'^-g(z')G(z; z') 



S B (x: z') + 



D-3 



Sc(x\ z') 



188) 



Interchanging x M and z^ gives, 



n 



rji'Ji 



6(x)A6(z)\ Q) = (A-type Terms 
-f d D x'yJ-g{x') G{x- x') 



S B (z:x') + 



D-3 



Sn(z:x') 



189) 



The term with two A#'s yields a surface integral of surface integrals that 
we can write as a volume integral of surface integrals using Green's 2nd 
identity, 

A9(x)A9(zf 1 



n 



n 



X 



d d 
S c (z; x)—G{x; x) - G(x; x)—S c (z; x) 



,(190) 



(D-3) 2 H 2 Jv 



d D x'J-g(x>) S c (z;x')n x/ G{x;x')-G{x;x')n x iS c {z;x') .(191) 



Of course we can use the identity J-g(x') n x >G(x; x') = id (x — x'), and 



the quantity J—g(x') \3 x >Sc(z; x') involves, 

-g{x') n x dA c (z';x') = id D (x'-z')+2(D-3)H 2 J-g{x') iA c {z'; x') . (192) 



The delta function in (j!92p gives another A-type Term whose form we work 
out in the next subsection. Hence we have, 



A0(z)A0(z)] fi) = (A-type Terms 



Sc(z\x) 2 



d D x'yj-g(x>) G(x; x')S c (z; x') .(193) 



(D-3) 2 H 2 D-3Jv 
The result is symmetric in x^ and z M so we can express it as, 

S c (x; z)+S c (z;x) 



n 



A6{x)A6{z) SI) = (A-type Terms 



+ f d D x'yJ-g{x') G(x; x') 



2(D-3) 2 H 2 

- S ^ X l + [j^f^G(z ] z')?^ . (194) 



D-3 
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Combining expressions 



with (11941) gives the desired form (I182|) . 



n 



A9{x)9(z)+9(x)A9(z)+A9(x)A9(z)]\q) = (A-type Terms 



S c {x; z)+S c (z;x 
2(D-3) 2 H 2 



d D x'yJ-g(x')G(x;x') 



S B (z: x') + 



S c {z; x') 



r v d D z'^/-g(z')G(z;z') 



S B (x; /)4 



S c (x; z' 



D-3 



D-3 
.(195) 



5.3 Absorbing the ^4-type Surface Terms 

We should begin this section by clarifying precisely what the A-type surface 
terms are. They reside entirely in the "Integral Term," and they consist of 
Sa/ (D — 3)H 2 plus the A-type surface terms induced by the gauge parameter 
A0L4]. We first reduce Sa to a pair of temporal surface terms, then derive 
similar expressions for the A-type surface terms from A0[A]. This will moti- 
vate our construction of the final gauge parameter 59 [A] which absorbs the 
A-type surface terms and gives full agreement with the invariant propagator. 
Recall that the surface integral Sa (x; z) is, 



S A (x;z) = d u -\c'J-g'g ,pa iAA l {x , - 1 z)& a G{x ] x')-G{x- 1 x')& a iA A '{x l - 1 z) , 

JdV J 

(196) 

where iAA>(x;z) is the derivative with respect to v (evaluated at v = (D — 
l)/2) of the Fourier mode sum, 



iA u (x; z) 



d D ~ l k 



jk-(x-Z) ] e(x -z°)u u {x , k)u* u (z°, k) 



+9(z°-x°)ul(x , k)u u (z°, k) .(197) 



Because G(x; z) is —i times the same mode sum (again evaluated at v = (D— 
l)/2) we see that the surface terms at spatial infinity make no contribution. 
One can therefore express Sa ( 

X) Zj clS db Fourier mode sum of temporal surface 

terms, 



S A (x; z) 



d D ~ l k 
(2u) D - 1 



x <^ u* A (x°, k)u* A (z°, k) x T{-kr] 2 ) - u A (x°, k)uA(z°, k) x ^(-krjx) L(198) 
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where rji and t] 2 are the initial and final times, respectively, and the function 
kri) is, 



7T f <9#W (-krj) dH^ (-kri) 



dvdrj 



(199) 



AHa 



dv 



drj 



H^^ftM) . (20 0) 



dvdrj 



D-l 
' 2 



This function ^(z) has the interesting property that it can be related to the 
product of two Hankel functions, without any derivatives with respect to the 
index or the argument [33J. To see the relation we define, 



S v {z) = z\HV{z) 



Quiz) 



Of course we have, 



d v H^\z)d z H^\z) - Hi 1 \z)d v d z Hj ) 1) 



7T 



H-kri) = -- x g v (-k V ) 



and the relation to E v is |33] 



2v 



d z Q v {z) = — r£ v (z) . 



(201) 
(202) 

(203) 
(204) 



The integration constant can be fixed using the asymptotic expansion for 
large z to give, 

poo f ( 

g v (z) = 2v dz'^-L. (205) 

J z Z 

The key identity for A0[A] to produce A- type surface terms is, 



d 2 y(x'; z') 
dx'°dz'° 



C'(y(x';z')) + 



dy(x'; z') dy(x'; z') 



dx'° 



dz'° 



C"(y{x'-z')) .(206) 



The A-type surface terms come exclusively from the delta function term; 
the other contributions produce 5-type and C-type surface terms we have 
already included. Note that because one gets a D-dimensional delta function, 
whereas the initial and final surface integrals are only (D — l)-dimensional, it 
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is necessary to regulate A#L4] to make the A- type surface term well-defined. 
An obvious regularization is to integrate the initial and final time surfaces 
over a small range of duration Arj = 2e, 



A9Jx) 



1 



2e(D-3)H 

,D-2 



dx'° - I dx 

rfi— e Jf]i—e 



xa°r 2 { — A (x')JLg(x;x') -G(x;x')^\— A (x')]\. (207) 
a r i ox' v ox' u l a T / J 



Let us now work out the A-type surface term from A9 t (x) x 9(z). The 
full expectation value is, 



n 



A6 e (x)9(z) 



dx'°- / dx'° 



jD-i-tt f jD j a x> 2 G(z\z') 
X Jv 2e{D-3)H 



d_ 

dz 70 



a°-HA c (x'; z')^ G(x; x>) - a^G(x; x')-JL %Ac (x'; z') ). (208) 



dx'° 



However, we already accounted for most of this in the previous subsection; 
it is only the delta function from using (12061) on the final surface term which 
makes the new contribution we seek, 



n 



A9 e {x)9(z) 



n 



—i 



dx'°- / dx'° 

772— 6 Jrji—e 



A-type 2e(D-3)H 
xfd^x'f d D z'G(z; z')a°- l G{x-x')5 D {x'-z') .(209) 



V 

The integration over z'^ is not affected by our regularization of A8 t (x) 

'■112 



d D z' = I dz' u \d u - v z" . 



/0 \jD-\ - 



(210) 



It is therefore only half the x'° range over which the delta function can be 
saturated. Taking the unregulated limit gives, 



e^0 



T*[Ae e (x)e(z)\ n 

= lim ~ l 

™2e{D-3)H 

i 



A-type 



V2 rvi+ e 
>dx'°- dx'° 

r) 2 -e Jrn 



d D -^G{z] xy^Gix; x') , (211) 



J\„D-1 



2(D-3)H 



d D -^ 



(212) 



a x , iAa(x; i')iAa(z; x') 
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One obviously gets the same result (12121) from 9{x) x A9 e (z) so the total 
for these "mixed" terms is, 



\im(nT* 



A9 e (x)9(z)+9(x)A9 e (z) 



n 



A-type 



(D-3)H 



a x , iA A (x; x')iA A (z; x') 



(213) 



s k < 3 ~*)> a°- l u\{x\ k)u* A (z°, k) u A (r, 2 , k) 



(D-3)H J (2tt) d - 1 



^^^^^u\(x\k)u\(z\k) x l[Hj»(-kn2) 



a°- l u A {x\ k)u A {z°, k) u A ( Vl , k)'\, (214) 



(D-3)H 2 



-u A {x°,k)u A (z°,k) x ^[H^i-k^Y^. (215) 



*2 



Expression (12 1 5[) combines nicely with the A-type surface term from 9, 

Sa(x;z) = i r d D - x k j£.{x-?) 
(D-3)H 2 (D-3)H 2 J (2n) D - 1 

x\u A (x ,k)u A (z°,k) x T(-kr} 2 ) - u A (x° , k)u A (z° , k) x T*(—kr]x) 1.(216) 



By partial integration we can express T(z) as, 



T(z) 



x 2v / dz'-, [hU(z') 

4 J z Z L 

71 



H^\z) 



2 7T r°° , 1 d 
~A x I T^dl' 



z'»HV{z>) 



So the first term of ( 12 18f) is cancelled by (12151) . 

The full expectation value for A9 e (x) x A9 e (z) is, 



tt T* A9 e {x)A9 e {z) 



/ dx'°- / dx'° 

'ri2—e Jr\\—t 



n 



-i 



Ae 2 (D-3) 2 H 2 
d D - x x' a D r 2 



dz'°- dz'° 

772 -e J-qi-e 



d D -^a D r 2 



(217) 
(218) 
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f-A / / /x dG(x; x') dG(z: z') ^. ..dG(x;x')diAc(x';z' 
x ^ iA c (x ; z ) — 7^ 5-^ G(z; z^ 



-G(x; x 



dx'° dz'° 
f .dG(z;z')diA c (x';z') 



dz'° 



dx'° 



+ G(x;x')G(z; z' 



dx'° dz'° 

rx dHA c (x>;z>) 



dx f0 dz'° 



•(219) 



As with the mixed term (I212p we have already reduced most of this in the 
previous subsection. The only new contribution derives from the delta func- 
tion one obtains by using (12061) on the final term, 



n 



A8 e {x)A6 e {z) 



n 



dx'°- / dx'° 



d D ~ x x' a D r 2 



x 



A-typc Ae 2 (D-3) 2 H 2 

d D - 1 ^ G(x; x')G(z; z')5 D (x'-z') , (220) 



Ae 2 (D-3) 2 H 2 

i 



4e 2 (D-3) 2 H 2 J (2tt) d - 1 



/ dz'° 


- / dz'° 




Jrii—e 




/ dx'° 


Jrj2—c 


Jr]l—€ 


r d°-\ 


ik-(x— 



ik-(x-z)) u *(r° lAn,*.(v Q 



(x\k)u\(z\k) d V a»- 2 u A (rj,k) 



-u A (x°,k)u A (z°,k) d V a°- 2 u* A (r],k) . (222) 



The integrations with respect to rj can be performed, but it is not possible 
to take the unregulated limit, 



V'2+C 



dr]a v u A {r},k) 
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dr)r)[HM(-kr]) 
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dz v { ) 



z=-k(r)2+e) 



(223) 

(224) 
, (225) 



z=— k(r)2— e) 



-- x 2e x m 
4 



HpH-krfr)] +0(e 3 ) . (226) 
Combining the various A-type surface terms gives a result of the form, 



(^4 — type Terms 



—i 



d D -\ 



ik-(x—z) 



{D-3)H 2 J {2ti) d - 1 
x \u A {x°, k)u* A {z°, k)A( V2 , k, e) - u A (x°, k)u A {z°, k)A*( Vl , k, e) ) .(227) 
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The function A(rj, k, e) is, 



71 



2e(D-3) 



+ dz' 



1 d 



+ 0(6) 



(228) 



We seek a homogeneous gauge parameter 89 [A] (x) which cancels ( )227H228p . 
in the limit that e goes to zero, without changing the "Other Term." If we 
construct it from Ai, rather than Aq then there will be no interference be- 
tween A0[A] and 59[A]. Suppose further that 50[A](x), like A9[A](x) involves 
an integral over a dummy variable x ,tJl , and that the field Ai(x') is differen- 
tiated with respect to x'°. What we want is that the expectation value of 
50 [A] (a;) x Ai(z) is zero unless there is also a derivative with respect to z°. 
If we can construct a 59 [A] with this property then the only nonzero contri- 
bution to the transformed propagator will come from 50[A](x) x 

It is simplest to construct the term we want by analogy with the simple 
harmonic oscillator, whose Heisenberg position operator is, 



\j2muj 



ae-^ + ate™ 



(229) 



Note that we can isolate the raising and lower operators by taking linear 
combinations of q and iuq, 



q(t)+iuq(t) 



2iu 



: a ] e iwt 



q(t)-iuq(t) 



—2icu 



a e 



—ICjt 



(230) 



\j2muj y/2mu 
We assume the usual commutation relations and ground state 

[a,a f ] = l , a\Q) = = (fi|a f . (231) 

If t comes before the last time t 2 , and after the earliest time t\, then we have, 



rp* 



(q{t 2 )+iuq{t 2 ))q{t) 
(q^-iuqih^qit) 








v t < t 2 , 

V tx < t . 



(232) 
(233) 



The only nonzero expectation value comes from the T*-ordered product of 
two factors of the combinations 



rjiif 



(q{t 2 )+iuq(t 2 ))(a{t' 2 )+iuq{t' 2 
q(ti) J (q(t[) -iuq^ 



ni 



5(t 2 —t' 2 
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(234) 
(235) 
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We construct a gauge parameter SO [A] (x) with the desired properties by 
analogy. The free field mode sum for A^x) is [TTJJ, ITS] , 



Ai 



d n ^k 



(x) = J ~ D "_ x <a x u B {x°, ky^frik) + a x u* B (x {> , k)e-^(3j{k)y (236) 



where the u B (rj, k) mode functions are given by f 1 1 1 8 j) with index v = (D — 
3)/2 and the f3i(k) are canonically normalized annihilation operators. One 
can isolate fli(k) by taking the spatial Fourier transform, 



Ai{x^k)= Id^xe-^Aix^x) 



(237) 



Now form linear combinations analogous to ( 12301) . 

1 ~ - d <9 r l~- 

- A{V, k )TrU* B (v, k) - u B (r), k)- Ai(rj, k) 

a or] or] l a 

1 ~ A d <9 rl ~ - 

- Ai(r), -k)—u B (r], k) - u B (r), k)— - Ai(r), -k) 
a or] or] l a 

It follows that the desired gauge parameter is, 
56 Ax) 



iPi{k) 
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(238) 
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ka x i 
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ka.ri 



1.(240) 



6 Discussion 

There are two generic ways to freeze local symmetries: 

• Exact Gauge Fixing, in which the fields are made to obey some equa- 
tion; and 
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• Average Gauge Fixing, in which a term is added to the Lagrangian. 

We have shown that certain average gauges cannot be derived from the canon- 
ical formalism on manifolds such as de Sitter for which there are linearization 
instabilities. Ignoring this problem in electrodynamics causes the vector po- 
tential to possess an unphysical and incorrect part which drops out of the field 
strength but affects interaction energies. This may be the origin of the on- 
shell singularities found in Feynman gauge for the one loop self-mass-squared 
of charged scalars on de Sitter [18J. 

We have also constructed the field-dependent gauge transformation that 
enforces exact, Lorentz gauge on de Sitter electrodynamics. This was ap- 
plied to the photon propagator from a non-de Sitter invariant, average gauge 
and the result agrees exactly with the de Sitter invariant solution previously 
obtained from solving the Lorentz gauge propagator equation [19J. It was 
already known from adding the compensating gauge transformation to the 
naive de Sitter transformation that the propagator in the non-invariant gauge 
shows no physical breaking of de Sitter invariance [10] . So the fact that our 
transformation technique produces an invariant result demonstrates that the 
technique indeed eliminates unphysical breaking of de Sitter invariance. 

In a subsequent work we will employ the same technique to transform 
the graviton propagator from a non-de Sitter invariant, average gauge [TP] 
to the exact and de Sitter invariant, de Donder gauge. Adding the com- 
pensating transformation shows that the breaking of de Sitter invariance in 
this propagator is physical [16], so the expectation is that the transforma- 
tion technique will not remove it. Because the transformed propagator will 
obey a de Sitter invariant gauge condition, this should settle the issue about 
whether or not free gravitons have any de Sitter invariant states. Note that 
simply obeying a de Sitter invariant propagator equation does not guarantee 
a de Sitter invariant solution, as the case of the massless, minimally coupled 
scalar proves [6]. Note also that physical graviton modes obey precisely the 
same equation as the massless, minimally coupled scalar [1] . 

Constructing the de Donder gauge propagator is a worthy goal in its 
own right for two reasons. First, exploiting the gauge condition makes a 
vast simplification in tensor algebra [15]. Second, using a de Sitter invariant 
gauge would preclude the need for noninvariant counterterms, even though 
the actual propagator is not de Sitter invariant [HI [T5] . 

A significant technical result of this paper is the "Convolution Identity" 
(11311) for integrating the propagator iA^ of a massless, minimally coupled 
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scalar up against the propagator %A V of a massless scalar with conformal 
coupling, 

t=r*\(5±f-S\. (241, 



D(D-l) 

The result follows from Green's second identity, 



iA u (x; z) —iAa(x; z 

a~ x \i —q[x' )il\a[x\x' )ii\ v {x \ z) = 

v 



-if d D x'^-g(x') iA A (x; x')iA u (x'; z) 



iA v (x'; z)d' a .iAji(x; x') —iAa{x\ x')d' a iA v {x' ] z) 



. (242) 



We expect this to be of great utility in the subsequent graviton project be- 
cause field dependent gauge transformations result in precisely such convo- 
lutions. 
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